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Transport through smaU interacting systems connected to noninteracting leads is stud- 
ied based on the Kubo formahsm using a Ehashberg theory of the analytic properties of 
the vertex part. The transmission probability, by which the conductance is expressed as 
g — (2e^//i) / de (— 9//9e) T(e), is introduced for interacting electrons. Here /(e) is the Fermi 
function, and the transmission probability T{e) is defined in terms of a current vertex or a 
three-point correlation function. We apply this formulation to a series of Anderson impurities 
of size N {— 1,2,3,4), and calculate T(e) using the order C/^ self-energy and current vertex 
which satisfy a generalized Ward identity. The results show that T(e) has much information 
about the excitation spectrum: T(e) has two broad peaks of the upper and lower Hubbard 
bands in addition to N resonant peaks which have direct correspondence with the noninteract- 
ing spectrum. The peak structures disappear at high temperatures. 
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§1. Introduction 

Effects of electron correlation on transport through small interacting systems have been a subject 
of current interest. For instance, the Kondo effectS) in quantum dots has been studied intensively, 
and recent experiments!' §11* have demonstrated a qualitative agreement with theoretical predic- 
tionsBi'B 

To investigate transport properties of small systems, it is essential to use theoretical approaches 
which are able to treat correctly both the interaction and quantum interference effects. The Keldysh 
formalismil' has been applied to nonequilibrium systems-BIll* Caroli et al have shown for nonin- 
teracting electrons that the current is determined by single-particle excitations between the two 
chemical potentials of the left and right leads.& However, to our knowledge, the same feature has 
been shown to hold for interacting electrons only in a special case where the hopping matrix el- 
ements, which connects an interacting region at the center and the two leads, satisfy a certain 
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condition. 0) Although this condition for the connectivity is not so important for a single Anderson 
impurity, it restricts the application of the formulation to interacting systems consisting of a num- 
ber of quantum levels. Our motivation of a series of works is based on a conjecture: as far as the 
linear response is concerned, the current through small interacting systems may be determined by 
low-lying states near the Fermi energy. It has been clarified without the assumption of the connec- 
tivity that the conductance is written in terms of the Green's function at the Fermi energy e = 0, 
at zero temperature T = 0. This is caused by a ground-state property of Fermi systems, i.e., 
the imaginary part of the self-energy vanishes at e = 0, T = 0.0^ Consequently, the contribution of 
the vertex corrections vanishes. The proof was provided based on the Kubo formalism making use 
of the diagrammatic analysis jil'lll* and it justifies an interpretation of the conductance in terms of 
the free quasi- particles0llll) 

The purpose of this paper is to generalize the formulation to finite temperatures. To this end, we 
have carried out the analytic continuation of the vertex part making use of a Eliashberg theory of the 
analytic properties, which was originally proposed to the translational invariant Fermi systemsEl^ 
and extended to lattice systems.El^ Compared to these examples, the situation we are considering is 
different in the following points: systems have no translational symmetry, and noninteracting leads 
connected to the interacting region keep the initial and final scattering states unrenormalized. One 
of the main results is that the conductance can be expressed as eq. (2.35) using the transmission 
probability T{e) which is defined in terms of a current vertex eq. (2.36) or a three-point correlation 
function eq. (3.6). It can be also written in the Lehmann representation as eq. (A-6). 

We apply the formulation to a series of Anderson impurities of size N (=1,2,3,4). For N > 2, the 
off-diagonal (non-local) elements of the self-energy and current vertex play important roles on the 
transport properties Therefore, it is essential to perform calculations with a current-conserving 
approach. We calculate T(e) in an electron-hole symmetric case using the order C/^ self-energy and 
current vertex which satisfy a generalized Ward identity. At low temperatures, the value of the 
conductance is determined by the low energy part of T(e) because of the Fermi factor —df/de in eq. 
(2.35). Nevertheless, the high energy part of T{e) also has much information about the excitation 
spectrum of interacting electrons: T(e) has two broad peaks of the upper and lower Hubbard 
bands in addition to N resonant peaks which have direct correspondence with the noninteracting 
spectrum. At low temperatures, the peak structures near the Fermi energy depend on whether 
is even or odd: the Kondo peak is situated at e = for odd while it is a valley for even A^. The 
even-odd oscillation becomes small with increasing T and disappears at high temperatures T > Tk, 
where Tk is the Kondo temperature. 

In §||, we describe the formulation, and introduce the transmission probability T(e) carrying out 
the analytic continuation of the vertex part. In the alternative expression of T(e) is given in 
terms of a three-point correlation function, and the current conservation is discussed in terms of 
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a generalized Ward identity. In we apply the formulation to a linear chain of the Anderson 
impurities. Summary is given in The Lehmann representation of T(e) is given in Appendix. 

§2. Conductance and Transmission Probability 

In this section we describe a finite temperature formulation for the current through small inter- 
acting systems connected to reservoirs based on the linear response theory. We perform the analytic 
continuation of the vertex part following Eliashbergjll' and introduce the transmission probability 
for interacting electrons, T(e), by which the conductance can be expressed in a Landauer type 
form.§ilii) 

2.1 Model 

We start with a system which consists of three regions; a finite central region (C) and two 
reservoirs on the left (L) and the right {R). The central region consists of N resonant levels, and 
the interaction Uj^j^-j^ji is switched on only for the electrons in this region. We assume that each 
of the reservoirs is infinitely large and has a continuous energy spectrum. The central region and 
the reservoirs are connected with the mixing matrix elements and f^, as illustrated in Fig. 1. 
The Hamiltonian is given by 

H = Hl + 'Hr + 'Hc + 'Hc'^ + Hmix , (2.1) 

= E E ( -tfi - M ) 4s-. ' (2-2) 

ijeL a 

^« = E E ( - '^'J ) 4c,-. , (2.3) 

'Hmix = ~ E ''^i ( 4crCocr + Cocr'^lcr ) ~ E ^-R ( ^N^Xa^Na + 4o-'^Af+lo- ) • (2-6) 
a a 

Here 4 ('^jo-) creates (destroys) an electron with spin a at site j, and ^ is the chemical potential, tfp 
tf-, and t'^j are the intra-region hopping matrix elements in each regions L, R, and C, respectively. 
The labels 1, 2, . . N are assigned to the sites in the central region. Specifically, the label 1 (N) 
is assigned to the site at the interface on the left (right), and the label {N + 1) is assigned to the 
site at the reservoir-side of the left (right) interface. We assume that the hopping matrix elements 
are real, and the interaction has the time-reversal symmetry: 114^3-21 is real and f743;2i = f^34;i2 = 
C^i2;34 = f42;3i = C^i3;24- We will be using units h = 1 unless otherwise noted. 
The single-particle Green's function is defined by 

Gjfiie^) = - rdr(r,c,.,(T)c],,(0)) e''-\ (2.7) 
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Here /3 = 1/T, = {'^'m + Cj^{T) = e'^^c^o-e "^^j and (• • •) denotes the thermal average 



Tr 



e-/3W . . . 



/Tre"^^. The spin index has been omitted from the left-hand side of eq. (E 
assuming the expectation value to be independent of whether spin is up or down. Since the 
interaction is switched on only for the electrons in the central region, the Dyson equation can be 
written as 

G,,{z) = G%{z) + GUz) G,,,{z) . (2.8) 

j'j'ec 

Here G^j{z) is the unperturbed Green's function corresponding to the noninteracting Hamiltonian 
Ti^ = Hl + Ti-R + TCl; + TCmix- The summations with respect to i' and / run over the sites in 
the central region, and is the self-energy due to the interaction TC'q*. Owing to the time- 

reversal symmetry of 7i, these functions are symmetric against the interchange of the site indices: 
Gij{z) = Gji{z) and T,ij{z) = Sjj(z). Note that at T = the single-particle excitation at the 
Fermi energy z = iO~^ does not decay, i.e., ImS^(O) = 0. We will treat z as a complex variable, 
and use the symbol + (— ) in the superscript as a label for the retarded (advanced) function: 
S±(e) = S,,(6±iO+). 

2.2 Analytic continuation of the vertex part 

We next consider the conductance based on the Kubo formalism. The conductance g is given by 
the w-linear imaginary part of a current-current correlation function: 



g = e hm , (2.9) 

a;-+0 lUJ 



2,. Kaa'{iO + iO+)-Kaa'{^0+) 

e lim : 

Kaa'im) = f^dr {T, Ja{r)Ja'm (2.10) 
Jo 

Here vi = 27r///? is the Matsubara frequency, and the retarded function can be obtained through 
the analytic continuation Kaa'{^ + '^0^) = -f^oa'(i^Olijy;-^c>j+io+- '^^^ current operator Jq, for q = L 
or i?, is defined by 

= i ^ -^ij ( cjv+i^CiVa - cLc7v+ia ) • (2-12) 

Here Jl is the current flowing into the sample from the left lead, and Jr is the current flowing out to 
the right lead from the sample. These currents and total charge in the sample, = J^jeCa 4o"So"' 
satisfy the equation of continuity dp^/dt + Jr — Jl = 0. Due to this property, the value of g given 
by eq. (2.9) is independent of a and q'.© Note that K^'ai^) = K^a'iz) because of the time-reversal 
symmetry of Ti. 

We calculate the w-linear part of Ko,a'{z) taking a and a' to be R and L. The function Kmi^i^) 
can be separated into the two terms as shown in Fig. 2: 

KrU^i.) = K^^liiu) + 4'i(i^) , (2.13) 
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/51 



i£ cr 



GN,i{is + ii^) Go,N+i{is) + Giv+i,o(i£ + ii') Gi,Ar(ie) 
-GAr+i,i(i£ + iv) Go,Ar(ie) - GAr,o(i£ + iz^) Gi,Ar+i(ie) 
(i^^L)(i^'ii) i EE E 

^ i£,ie' crcr' {i}eC 

(ie, ie'; iz^) 



(2.14) 



(2.15) 



Here FJ^'j^. ^^^^(ie, ie'; ii^) is the total vertex part due to , and the summation with respect to 
{j} runs over the sites in the central region Ji, • • • ,i4 G C. Since the leads are assumed to be 
noninteracting, the Green's function satisfies the following relations at the two interfaces: 

Goj(z) =-gL{z)vLGij{z) for 1 < j < iV + 1 

Gj^N+i{z) = -Gj,N{z) gjiiz) for < j < iV 

Here g^iz) and gji{z) are the local Green's functions at the interfaces of the isolated leads, 
and N + 1, respectively. Using these relations, KjiL{ii/) can be rewritten in terms of the Green's 
functions in the central region as 



(2.16) 



KRL{iu) = 4EE E AL(ie,i£ + iz^)Gij4(ie)AR;j4jj(i£,ie + iz/)Gjii(ie + ii/) 

^ 16 cr jij4eC 

Here the bare and full current vertices are defined by 

AL(i£,i£ + iiy) = -ivl [g^(i£ + ii/) - gi(i£) ] , 
AR(i£, i£ + w) = ivji [gfl(i£ + iiy) - g^iie) ] , 

^i4.i(i£,i£ + i^)= ^EE E rj,l;,3,,(i£,i£';iz^) 

ie' a' iaiseC 

X Gj^Ni^s') AR(i£', ie' + ii/) G^j^ i}^' + i^) ■ 



(2.17) 

(2.18) 
(2.19) 
(2.20) 

(2.21) 



The spin index has been omitted from the left-hand side of eq. (2.21) since it is independent of 
cr due to the rotational symmetry of the spin space. The w-linear part of the retarded function 
KRL{ijo + zO"^) can be obtained from eq. (2.17) by carrying out the analytic continuation: 

KRL{u + iQ+) = Y, E 
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-oo 27ri 



f{e + cv)-f{e) 



(6, e + ^) Gy, (6) A^.^.^ (6, 6 + ^) G+ 1 (e + ^) 



J2] 



+ /(^ + ^) 4"(e,e + ^)Gr,,(6)Ag|.^^.^(6,e + a;)GTi(e + ^)|. (2.22) 

Here /(e) = [e^*^ + 1]^^. The superscript [k] for A; = 1,2,3 is introduced for the current vertices in 
order to specify the three analytic regions of Aji-j^jj^{z, z + w) in the complex z-plane, which are 
separated by the two lines, Im (z) = and Im {z + w) = 0, as shown in Fig 3. The explicit form of 
the bare current vertex is given by 



for a 



L, R. Thus, in the limit of uj 



0, AL''(e,e) 



for k 
for k 
for k 



^L"(e,e) 



1 

2 
3 



(2.23) 



0. Furthermore, Xj^ {e,e) 



and Ao 

2rR(e) and X^l'{e,e) = —2TL{e) with Ta{^) = — Im [wq^oI^)]- order to calculate the w-linear 

\k] 

part of eq. (2.22), we next consider the behavior of Pj^j^(e,e + lu) for small lo carrying out the 
analytic continuation of Pj_^j-^{ie,ie + iv) using eq. (2.21). It can be done following the Eliashberg 
theory of the analytic properties of the vertex part© As a function of the complex variables, 
Tj^j^.j^j^{z,z';w) has singularities along the lines shown in Fig. 4: the vertical lines Im (z) = 
and Im {z + w) = 0, the horizontal lines Im (z') = and Im (z' + w) =0, and the diagonal lines 
Im (z — z') = and Im (z + z' + w) = 0. These lines divide the plane into 16 regions, each of 
which corresponds to a function of r'j^j^.j,^j^{z,z';w) which is analytic in that region in any of its 
arguments. The rectangular regions in Fig. 4 are numbered by indices [k,k'] for k,k' = 1,2,3. 
Some of these regions are divided into parts by the diagonal lines, and these parts are denoted by 
Roman numbers. We will use this notation also for the vertex part in eqs. (2.25)-(2.33). Rewriting 
the summation over the Matsubara frequency ie' in eq.(2.21) in terms of a contour integral along 
these lines and then performing the analytic continuation, Pj'^l-^i^, e + lo) is obtained as 




+ 



de' 



[3], 



UJ 



. (2.24) 
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Here £^'^j.'^'';](e, e'; w) is defined by 



Jlj2;j3j4^ 

paa' [11] / / \ 



= -[/(e' + a;)-/(eO]r-5iS..(^>^';-i. 

= + e + 0.) [ r^SS -) - r^:S? -) 



^o-cr' [12] / / 



/^itit' [13] / / 



^o-cr' [21] I I 



aa' [22] 



[23] / / 



JIJ2; J3J4 
^aa' [31] / / 
^ili2;i3j4i^'^ 



cro-' [32] 



J1J2;J3J4^^' ^ ''^ 

^itit' [33] / / 



'[12] 



-.ctct'[22:IV1 



, T>U ' \ rT^<7<7'[22;II]/ / X ^(t(t'[22;III] 

+Vb{e' + e + a;) 



/(^')r-2S"'(^'^';-)-/(^' + -)r;:^2;.3.4^ 

(e,e';a;) 

paa'[22;III]. , x paa'[22;IV]. , x 

= + e + [ r-^S];i ie,e';.)- r-^Sll ^'^ 

= -[/(e'+.)-/(.o]rj-;i:5^^(.,e';..), 



(2.25) 
(2.26) 

(2.27) 

(2.28) 



(2.29) 
(2.30) 

(2.31) 
(2.32) 

-Jij2;j3j4v-'- • (^•^^) 
Here V denotes the Cauchy principal value, and 6(e) = [e^*^ — The Fermi and Bose distribution 
functions came through the summation over ie'. We can now examine the behavior of pj^j^(e, e + 
a;) for small lo. The first and third terms of the right-hand side of eq. (2.24) vanish at w = 
because A^^(e,e) = and A^'(e, e) = as already mentioned. Furthermore, due to the factor 
[f{e' + u)-fie')] appearing in eqs. (2.26) and (2.32), £f;^^^J^^^(e,e';0) = /:^^^J^j^(e,e';0) = 
0. Therefore in the Hmit of w ^ 0, Pjl\,{e, e) = 0, Pj%ie, e) = 0, and 

= 2E E <S4(^'^';0) Gj,M^')^n(^')G%,^e') . (2.34) 

Similarly, through eq. (2.20) we have, A%,j,{e,e) = 0, A%,j,{e,e) = 0, and A%,j,{e,e) = 
2Tii(€) dNji^Nji + Pj4,jii^'^)- Therefore, the a;-linear part of Krl{^ + ^0"'") comes only from the 
second term of the right-hand side of eq. (2.22), and thus the conductance can be expressed as 



2e^ 
~h 



de 



9/\ 
de 



,[2] 



r(6)= E 2r^(e)Gr,,(6)A^;.^.^(6,6)G+,(6). 



(2.35) 
(2.36) 



Among nine functions of C-^j'^j.l^j,^]^ given by eqs. (2.25)-(2.33), only one function ^J^jj- j^j4 I'elevant 

[21 f [21 1 * 

to the conductance through eq. (2.34). We note that A^^^-,(e,e) = | A|^^,j (e, e) | due to the time- 
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reversal symmetry (see also Appendix), and thus 7"(e) is real. The contributions of the bubble and 
vertex diagrams can be separated as T(e) = T^^^e) + T^^\e), 

T(-)(e) = 4TUe)G^^{e)TR{e)G+,{e) , (2.37) 
tW(6)= 2ri(6)Gy^(6)P,[3^(e,6)G+,(6). (2.38) 

At low temperatures, the conductance is determined by the behavior of T(e) at small e because 
of the Fermi factor —df/de in eq. (2.35). Specifically, the vertex contribution vanishes at zero 
temperature.© This can be reexamined as follows: at w = 0, e = and T = 0, the [2,2] vertex 
function vanishes, i.e., ^J^jj; j^i4(^' ^) ~ ^' ^° restriction of the phase space by the 
distributions /(e') and 6(e') in eq. (2.29). Thus from eqs. (2.34) and (2.38) we have, pjj^(0,0) = 
and T('')(0) = 0, at zero temperature. 

In this section, we have assumed the time-reversal symmetry of Uj^j^-j^j^, but not assumed the 
precise form of this interaction. Thus, the present formulation can be applicable to the interaction 
which vanishes smoothly at the interfaces. We have also assumed that the two leads are nonin- 
teracting. Due to this assumption, we have no wavefunction renormalization in eq. (2.36). If the 
interaction is switched on also in the leads, the wavefunction renormalization may be necessary to 
define the initial and final scattering states .111* 

§3. Current Conservation and Ward Identity 

In this section, we provide an alternative expression of the transmission probability T(e) in terms 
of a three-point correlation function, i.e., eq. (3.6). This expression may be suitable for nonper- 
turbative calculations such as numerical renormalization groupii'0@^ and quantum Monte Carlo 
methodB The Lehmann representation of T(e) can be obtained using the correlation function, 
and it is given in Appendix. We also discuss the current conservation in terms of the generalized 
Ward identity eq. (3.7). This identity can be rewritten in the form of eq. (3.12), and it shows a 
relationship between the self-energy and current vertex. 

The three-point correlation functions of the charge and currents are defined by 



^C;,y(T;Ti,T2) = (t,5pc(t)c^-,(ti)4,(t2)) , (3.1) 

$L;ii'(r; Ti, T2) = (r, Jl{t) C^Jn) c],^{t2)) , (3.2) 
<J>fi;i/(T;Tl,T2) = (r, Jk(t)c,,(ti)4,(t2)) , (3.3) 



where 6p(j = — (pq). The Fourier transform of these functions are given by 

<^>r.y n, ra) = ^ E '^rM^^^ + i^) e--(^^-^) e'^ , (3.4) 

ie,iu 

for 7 = C,L, R. In what follows, we will discuss the correlation functions in the central region 
assuming jj' £ C. The right-hand side of eqs. (3.2) and (3.3) contain the operators with respect to 
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sites in the reservoirs, and + 1, through Ji and J/j. The contributions of these two sites can 
be rewritten in terms of the Green's functions with respect to the adjacent sites 1 and N using eq. 
(2.16) with the Feynman diagrams for <^^.jj/(ie, ie + iz^). Therefore, these correlation functions can 
be expressed in terms of the current vertex introduced in §0 as 



i4jiec 



(3.5) 



Here Afi-j^j^ is the current vertex for J/j given by eqs. (2.20) and (2.21). The functions A^j^j-^ and 
■jiji are defined by the similar expressions; the bare vertices are given by AL(ie,ie + iz^) 5ij^ 5ij-^ 
and (^j4ji, respectively. Thus, the transmission probability, eq. (2.36), corresponds to the analytic 
continuation of eq. (3.5) in the region, ie + w ^ e + uj + iO''" and ie ^ e — zO"^, as 



Tie) 



(3.6) 



[21 

The Lehmann representation of $^l^j^(e,e) is presented in Appendix. 

In the rest of this section, we discuss the current conservation in terms of these correlation 
functions. Using the equation of continuity in the Matsubara form — (d/dr) 5p(j + i J/j — i Jl = 0, 
the generalized Ward identity is obtained aJ 



) 



d 



~ q:^ ^C;jj' {t; n,T2) + i ^R-jf (r; Ti , r2) - i ^i-jf {t; n , r2) 



(3.7) 



= (5(r - r2 ) Gjjf (n , r) - 6{ti - r) Gjj> (r, T2 ) . 

The Fourier transform of this identity can be expressed, using N x N matrices, as 

ii/#(7(i£; ie + ii^) + i^ij(i£; ie + iz^) — i^L(ie, ie + ii^) = G{ie) - G{ie + w). (3.8) 

Here G{z) = {Gjj'{z)}, ^^{z,z + w) = {^j-jj'{z, z + w)}, and A^{z,z + w) = {A^/-jj'{z,z + w)} 
for jj' E C. In this notation, eq. (3.5) is written as ^^{z, z + w) = G{z) A^{z, z + w) G{z + w). 
Therefore eq. (3.8) can be rewritten as 

iiy Ac{ie, ie + w) + i AR{ie, ie + iz^) - i ^^(ie, ie + iz^) = {G(ie + iz^)}""^ - {G{ie)y^ . (3.9) 



Furthermore, the Dyson equation eq. ( |2.8D can be expressed as {G{z)} ^ 
E{z): 



z 1 



c 



^21 



"^22 ~ ^ 



(3.10) 




















v\gR{z) 



(3.11) 
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and X'(z) = The mixing matrix Vmix has only two non-zero elements, because among the 

sites in the central region only the two sites at 1 and are connected to the reservoirs. Note that 
eq. (3.9) shows a relationship among the self-energy and current vertices, which has to be satisfied 
in approximate calculations in order to obtain current-conserving results. We now consider the 
analytic continuation of eq. (3.9) in the region ie + ii/ ^ e + a; + iO^ and ie — > e — iO"^. In the limit 
of a; —5- 0, it is given by A^^{e,e) — A^l\e,e) = — 2ImV,^j^(e) — 21m. Z!~^{e). This relation can be 

[21 

rewritten further, using eq. (2.34) and the corresponding expression for A^^J^.^^.^(e,e), as 

E E r ^^JiSU^'^^o)[G73^(^OrM60G+^,(60 + G-3,(eO 

We note that eq. (3.12) can be regarded as a relation between the quasi-particles damping and 
transport relaxation time. At T = 0, the damping rate due to the inelastic scattering vanishes 
ImEj^^-^ (0) = 0, and it is connecting to the property of the vertex part ^'jijj^jjj^i^^ ^'s 0) = 0- 
Specifically, in the single impurity case = 1, eq. (3.8) is a scalar equation with respect to 

[21 f2l 

the impurity site. It is written, performing the analytic continuation, as (e, e) — {e,e) = 
G~(e) — G~^{e). Furthermore, if the mixing terms satisfy the condition Tile) = XTfi{e), we have 
$P](e,e) = -A$g'(e,e). Then the conduct ance can be written in a well-known expression; 

..... = ^ ly l-l^G^U)] . (3.13) 

§4. Application to the Anderson Impurities 

In this section we apply the formulation described in previous sections to a linear chain of the 
Anderson impurities. A schematic picture of this model is shown in Fig. 5: the system consists 
of a series of impurities at the center and two reservoirs. This system may be considered as a 
model for a series of quantum dots or atomic wires of nanometer size. For = 1 and N = 2, the 
conductance has already been calculated with efficient numerical methods such as the numerical 
renormalization groupH) and quantum Monte Carlo methods. Our interest here is how the 
transmission probability T(e) depends on the frequency and temperature. In the previous work 
we have calculated the conductance using the order [/^ self-energy at T = 0, and found that 
the off-diagonal (non-local) elements of the self-energy play an important role on the transport 
properties .00-* In order to calculate the conductance at finite temperatures, the contribution of 
the vertex part eq. (2.38) has to be taking into account carefully. In the present work, we calculate 
all the matrix elements of the order U"^ self-energy and current vertex which satisfy the generalized 
Ward identity eq. (3.12) to obtain T{e) for A^ = 1,2,3, and 4. We note that the second-order 
perturbation theory has been applied to transport through the single Anderson impurity0'@S0* 
and related systems0'ii) in somewhat different contexts. 
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In this section the parameters of the Hamiltonian defined by eq. (2.1) are given exphcitly as 
follows. For tfj, we take the nearest-neighbor elements to be t and all other off-diagonal elements 



to be zero. We assume U. 



to be an onsite repulsion U. Furthermore, we will concentrate on 
the electron-hole symmetric case taking the parameters to be /U = and + U /2 = where eo is 
the onsite energy — t^' = eo. Then the Dyson equation can be written in the N x N matrix form 
as, {Giz)}-^ = {G-{z)}-^-Uiz), 

e-vlgl{e) t 
t e t 



{G°(e + iO+)} ' = 



(4.1) 



t e t 
t e-vlgl{e) 

Here v1^g^{e) for a = L,R corresponds to the mixing term eq. (3.11) due to the coupling with 
the reservoirs, and the imaginary part Ta{e) = — Im g+(e)] is finite because each reservoirs 
has a continuous energy spectrum. Thus, the unperturbed Green's function G°{z) describes a 
system of N resonant states which have finite level width. Note that the Hartree-Fock term of U is 
already included in eq. (4.1), and thus ^{z) is the self-energy due to the interaction Hamiltonian, 



+ ni|)/2], where rii, 



We now calculate the self-energy and current vertex in the second-order perturbation with respect 
to Hq"^. The order U"^ self-energy is shown in Fig. 6, and is given by 

'^7X2 



ImS;,,(e)=(- 



dei / de2 ImG°+(ei)ImG0+(e2)ImG°+(e2 - ei + e) 



3 3' 



33 



X /(62) [ 1 - /(ei) ] [ 1 - /(e2 - ei + e) ] + [ 1 - /(e2) ] /(ei) /(e2 - ei + e) , (4.2) 



|ReG;+(ei)ImG°+(e2)ImG°+(e2-ei + e) [1 - /(e2 - ei + e) ] /(e2) 
- ImG;+(ei) ImG°+(e2) Re G°+(e2 - ei + e) /(ei) /(e2) 



3 3 ' 



33 



■ ImG;+(ei) Re G°+(e2) Im G°+(e2 - ei + e) /(ei) /(es - ei + e) 



3 3 ' 



33 



(4.3) 



The order current vertex Pj^j^ (ie, ie-l-ii^), which satisfies the generalized Ward identity eq. (3.9) 
with this order U'^ self-energy, is shown in Fig. 7. The contribution of the four-point vertex part is 
given by 



Sjiji ^jiji U ^3431 (i^ - i£ ) 

a-' Sj,jo Sj-ijA U Xj.j,(ie — ie 



+ 6a-a' Sj^js 6j2j4, Yj^j^ (ie + ie' + iv) , 



(4.4) 
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(i^) = - 75 E Gin + i^") (i^") > 
(i-) = - ^ E ^Li (- - i^") ^L. (i^") ' 



(4.5) 
(4.6) 



where each terms in the right-hand side of eq. (4.4) corresponds to each diagrams in Fig. 7. The 
analytic continuation of Pj^j^ for ie + iz/ — > e + a; + iO"'" and ie ^ e — iO"'" can be performed exphcitly 
using this four-point vertex. That yields the expression corresponding to eq. (2.34) with 

-b{e'-e)\xt.{e-e')-XT (e-e') 



X I /(^') Kn (e + e' + u)- f{e' + uj) Yf. {e + e' + u; 



+ 6(e + e' + cu) [v+^^ie + e' + co) - Yr^^^{e + e' + 
Here X'^^-^{uj) and Y^j_^{uj) is the analytic continuation of eqs. (4.5) and (4.6), respectively, 

/OO 
- [ /(6) - /(e + a;) ] Im G^, {e + u) Im Gjt (e) , 
-OO 

/OO Af: 
- [ /(e) - /(6 - o;) ] Im Gjt (,, _ e) Im Gjt (e) . 
-OO 'I 



(4.7) 



(4.8) 
(4.9) 



Using eqs. (4.7)-(4.9), the analytic continuation pj^j^ can be obtained up to order U'^ as 



V vr / J —OO J —OO 
xImG?t (,2) ImG;+^(.2-.i+6) 

X |/(e2) [l-/(ei)] [l-/(e2-ei + e)] + [l-/(e2)]/(ei)/(e2-ei + e)| 

(TT\ 2 (.00 poo 
- / dei/ d62G;;^(e2)G°+.^(e2)r,,(62) 
/ J— OO J —OO 

XI-G;+ (^l) ImGjt (62-61+6) 

X |/(e2) [l-/(6i)] [l-/(e2-6i + e)] + [l-/(62)] /(ei) /(62 - 6i + e) | .(4.10) 

We note that the first and second diagrams in Fig. 7 give the same contributions. Further- 
more, the Green's functions have some special properties in the electron-hole symmetric case: 
Gtj{-e) = (-l)-^+^G-.(6), S±(-6) = (-l)-^+iS-.(6), and P^\-e,-e) = {-If-i pf{e,e). 
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Then the contributions of the second and third diagrams cancel each other out, and eq. (4.10) 
simphfies as 

Pf},S^,e) = 2 f^)' r dei r d62 G°;^(6i)G°4(6i)r«(6i) 
xlmGjt (£2) ImG°;,(e2-ei+6) 

X |/(e2) [l-/(ei)] [l-/(e2-6i + e)] + [l-/(62)] /(ei) /(62 - 6i + e) |.(4.11) 

Note that the restriction of the integration region due to the Fermi functions has the same form 
with that in eq. (4.2). Therefore, pj^|-^(e, e) and ImE^.,(e) show the similar T and e dependences. 

We calculate the transmission probability T(e) as follows: for the pair of the full Green's functions 
Gfj{e) appearing in eqs. (2.37) and (2.38) we use the interacting ones obtained through the Dyson 
equation with the order U"^ self-energy, and for the current vertex pj^j^(e, e) in eq. (2.38) we use 
the order C/^ one, i.e., eq. (4.11). The generalized Ward identity eq. (3.12) is satisfied in this 
treatment. We calculate all the matrix elements of the order U"^ self-energy and current vertex 
numerically using the unperturbed Green's function G^p{e) obtained through eq. (4.1). For the 
reservoirs, we assume that the local density of states at the interfaces to be flat and the bandwidth 
to be infinity. Then ^■^(e) is pure imaginary and independent of the frequency: ^^^■^(e) = — iFo 
for a = L,R. Furthermore, we concentrate on the case = (= F), where the system has an 
inversion symmetry. The results of the transmission probability T(e) for N = 1,2,3, and 4 are 
plotted in Figs. 8-11 for three values of U, where (a)~(d) correspond to the results at four different 
temperatures. In the figures for = 1, the frequency e is measured in units of A (= F/, + F/j) 
which is now equal to 2F. For > 2, we measure the frequency e in units of t, and take the mixing 
parameter to be T/t = 0.75. 

In Fig. 8, the transmission probability for the single Anderson impurity is shown for several values 
of [//(vrA): (— ) 0.0, (-o-) 1.0, and (-•-) 2.0. The temperature T/A is taken to be (a) 0.0, (b) 0.1, 
(c) 0.33 • • •, and (d) 1.0. In this case the transmission probability is written in terms of the Green's 
function: T(e) = — AlmG^(e) as deduced from eq. ( |3.13| ).ll3ij^ Therefore, Fig. 8 corresponds to 
the results of Yamada&YosidaB and Horvatic, Sokcevic & Zlatic3 the Kondo resonance at e = 
becomes sharp with increasing U, and two broad peaks which have an atomic character appear at 
±[7/2 for large [//(vrA) ^ 2. The Kondo peak decreases with increasing T and disappears at high 
temperatures in the cases of large U. We note that the spectral function of the single Anderson 
impurity has been calculated accurately with the numerical renormalization group method.il-* We 
have provided the perturbation results just for comparisons. 

In Figs. 9-11, the transmission probability for A^ = 2,3,4 are shown for three values of U/{2TTt); 
{—) 0.0, (-0-) 0.5, and (-•-) 1.0. The temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.33- • ■, 
and (d) 0.5. At low temperatures T(e) for each systems has A^ resonant peaks, which have direct 
correspondence with those of the unperturbed system. In addition to these A^ resonant states. 
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two broad peaks appears at ±[7/2 in the cases of large U. At T = the resonant peaks become 
sharp with increasing U, and valleys become deep. However, the height of the peaks decreases with 
increasing U. One exception is the Kondo resonance situated at the Fermi energy for = 3. The 
height of the Kondo peak is unity for any values of U, and it causes a perfect transmission. This 
is a general property of the odd N systems, and occurs if the systems have the inversion symmetry 
Tl = Tji together with the electron-hole symmetry.0) The characteristic energy scale is determined 
by the width of the Kondo peak Tx, and it decreases with increasing N. On the other hand, the 
transmission probability of the even N systems shows a minimum at e = 0, at low temperatures. 
The characteristic energy scale is determined by the width of the valley which eventually tends to 
the Mott-Hubbard gap in the limit of large N. Comparing the results for A > 2, we see that the 
high energy part at |e| ^ 2t shows the similar slope in the case of U/ (27rt) = 0.5. For this values of 
U, the levels of the atomic character ±17/2 stay inside the one-dimensional band edge ±2t which is 
well-defined for large A. The two atomic levels go outside of the edge in the case of U/ (27rt) = 1.0, 
and the high energy part of T(e) for A^ = 3 and A^ = 4 show the similar behaviors at |e| ^ 4t while 
that for A = 2 shows somewhat larger values. Therefore, as far as the high-energy behaviors are 
concerned, the A = 3 and A = 4 systems seem to capture some aspects of the large A^ systems. 
The low energy part of T(e) is rather sensitive to the temperature. In Fig. 9-11, we see that the 
transmission probability for low-energy excitations at —2t ^ e ^ 2t decreases with increasing T, 
and the structures of the resonance peaks washed away at high temperatures. For instance, at 
T/t = 0.5 the transmission probability for the A^ = 3 and A^ = 4 systems behaves almost the same 
in the whole range of e in the case of U/{2iit) = 1.0. Thus, at high temperatures, the even-odd 
oscillation disappears. 

The transmission probability T(e) can be separated into the bubble and vertex contributions, 
i.e., T^°'\e) and T(^)(e) defined by eqs. (2.37) and (2.38). We next show how each of these two 
parts contributes to the transmission probabilities discussed above. In Figs. 12-15, the bubble (o) 
and vertex (•) contributions are plotted for A^ = 1,2,3, and 4. In these figures, the transmission 
probability T(e) is also shown (solid line), and (a)-(d) corresponds to the results at different 
temperatures. The value of U is taken to be [//(vrA) = 2.0 for A = 1, and U/{2Trt) = 1.0 for 
A > 2. In the case of the single Anderson impurity, T(e) can be expressed in terms of the Green's 
function G^{uj) as mentioned. We have done the separation also in this case for the purpose of 
comparisons. As discussed in §p, the vertex contribution T^^^e) vanishes at e = 0, T = O.^-* In the 
present systems, T^^\e) shows the quadratic energy dependence for small e, i.e., T^^\e) oc at 
r = 0. This is a property of the Fermi-liquid and relating to the low-energy behavior of the quasi- 
particle damping, i.e., — ImS^j^(e) oc e^. The low-energy behavior of T(e) is mainly determined 
by the babble contribution at low temperatures. However, T(")(e) decreases with increasing T. For 
instance, in the A^ = 3 and A = 4 systems the bubble contribution almost vanishes at T/t = 0.5 in 
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the whole range of e, and the vertex contribution determines the total transmission probability. In 
the single impurity case, however, the bubble contribution remains finite even at high temperatures. 
This difference is caused by the fact that the bubble contribution T(")(e) = 4r^|G'^]^(e)p is a local 
quantity for = 1 while it is an inter-site correlation for N > 2. In the N > 2 systems, the vertex 
contribution T('')(e) dominates the high energy part of the transmission probability, at |e| ^ 2t, 
even at low temperatures. The low energy part of T^^\e) increases with the temperature. The 
value at e = shows the dependence at low temperatures, and it is also connecting with the 
temperature dependence of the damping rate ImS^^-^(O) through eq. (3.12). 

We calculate the conductance using the results of T(e). At low temperatures the conductance 
is determined by the low energy part of T(e) because of the Fermi factor —df/de in eq. (2.35), 
although the high energy part also has important information about the excitation spectrum as 
discussed in the above. In Fig. 16 the conductance is plotted as a function of the temperature 
T/t. Here the onsite interaction U/{2'Kt) is taken to be (a) 0.5 and (b) 1.0 for the systems of 
N = 2 (□), = 3 (•), and = 4 (O). The noninteracting results are also plotted for N = 2 
(dotted line), = 3 (solid line), and A^ = 4 (dashed line). In the case of A^ = 3, the perfect 
transmission occurs at T = 0, and the conductance decreases with increasing T. This is caused 
by the Kondo resonance, and is a common feature of the odd A^ systems. Since the Kondo peak 
becomes sharp with increasing U, the temperature dependence becomes steep for large U. We note 
that the conductance through the single impurity, which has been obtained with the numerical 
renormalization group, §§§) shows the similar temperature dependence. On the other hand, the 
conductance for even A^ shows a minimum at T = because of the valley structure of T(e) around 
the Fermi energy. Due to the reduction of the transmission probability at low energy part seen in 
Figs. 9 and 11, the conductance decreases with increasing U. Comparing the results for A^ = 3 and 
A'^ = 4, we see that the conductance shows almost the same T dependence at high temperatures. 
Especially in the case of U/{2iTt) = 1.0, the two curves almost overlap each other at T/t ^ 0.5. 
This is because the peak structures of T(e) are washed away at high temperatures as seen in Figs. 
10 and 11. This example seems to show an essential feature of the disappearance of the even-odd A^ 
dependence. The crossover temperature may be determined by Tk, and it decreases with increasing 
A^. The even-odd oscillation of the conductancJll* occurs at low temperatures T ^ Tk- 

§5. Summary 

In summary, based on the Kubo formalism we have introduced a transmission probability T(e) 
for interacting electrons connected to noninteracting leads, which is given by eq. (2.36) and derived 
using a Eliashberg theory of the analytic continuation of the vertex part. Among the 16 analytic 
regions of the vertex part Tj^j^.j^j^{ie,ie';w), the transmission probability can be expressed in 
terms of the analytic function of the [2,2] region - (e, e';Lij). It is defined by eq. (2.29) and 

obtained through the analytic continuation for ie -|- ii^ — s- e + ijO + iO"'' and ie — > e — iO"*" . Alternatively, 
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T(e) can be expressed in terms of the three point correlation function of the current ^^.jj,{e, e + m) 
as eq. (3.6), and it is also written in the Lehmann representation eq. (A-6). 

We apply this formulation to a linear chain of the Anderson impurities of size (= 1, 2, 3, 4) in the 
electron- hole symmetric case. We calculate T{e) using the order U"^ self-energy T,jj/{io) and current 

[21 

vertex Pjj,{e,e) which satisfy the generalized Ward identity eq. (3.12). The current conservation 
is fulfilled in this approximation. The off-diagonal (inter-site) elements of these functions play an 
important role on the transmission probability. At low temperatures T(e) has A^ resonant peaks, 
which have direct correspondence with the spectrum of the unperturbed system. For large U, T{e) 
has two additional broad peaks at zbC//2. The resonant peaks become sharp with increasing U . 
The peak structures are washed away at high temperatures, and the difference between the even 
and odd A^ systems becomes invisible. 

The formulation described in this paper can be generalized to multi-channel leads following along 
the similar line. The extension to the nonlinear response is also interesting, and it has been done 
for an out of equilibrium Anderson model up to the third order with respect to a bias voltage, ii) 
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Appendix: Lehmann Representation 

[21 

In this appendix we give a Lehmann representation of <I>^yj'(ie, ie + iv) and $^.^.^.,(e,e + uj). 
Specifically, we will concentrate on the diagonal term j = j' , and start with the Fourier transform; 

J^^ dr dn dra e'^'e^^^^e-'^'^^ (t, J(r) c(ri) c\t2)) = (3 5e+u,e' ^(ie, ie + ii^) ■ (A-1) 

Here we have suppressed the subscripts to simplify the notation. The integrations can be done 
explicitly inserting a complete set of the eigenstates, H\n) = En\n), aslil) 

$(ie, ie + iz^) = {l\c^ \m) {m\J\n) {n\c\l) 

Imn 

(ie + iu + E^- Ei){iu + S„ - En) 



~ (ie + En- Ei){ie + iu + E„,- Ei) 

~ (ii/ + Em- En){ie + En- El) 



16 



+ \Y.'^Mn){n\J\m){m\c^\l) 



Imn 



-I3E„ 



{ie + Ei-Er,){iv + Er,-Er, 



+ 77 



(k + Ei- Er,)(ie + {u + Ei- E^) 



(A-2) 



(k + 11' + El- Em){iiy + En- Em] 
Here Z = Tie~^^. We now introduce the spectral functions defined by 

A^{loi,U2) = e"'^^' {l\c^m){m\J\n){n\c\l) 6{iui + El - Em) S{U2 + El - En) , (A-3) 



Irnn 



A^{uJi,U2) = -^Y. e"''^' {l\c\n){n\J\m){m\c^\l) 5{uJi + Ei - E^) 5{u2 + Ei - E^) . (A-4) 

Imn 

Using these spectral functions, eq. (A-2) can be rewritten as 

$(i£, ie + iz^) = - / da;ida;2 ^-(wi, 0^2) 
^ J— 00 

g-/3a;i _ Q-pu>2 e~^'^i + e~^'^'^ + 2 



+ 



+ 



ie + CJ2 ie + iu + cJi J iv + uji — L02 (ie + 0^2) (is + 11^ + ^01) 
1 

- / du;ida;2^+(wi,W2) 
J— 00 



1 



+ 



ie — L02 ie + 'w — ojiJ iu + UJ2— (ie — ^2) (ie + 'w — ui) 



■ (A-5) 



Carrying out the analytic continuation ie — >■ e — iO"*" and ie + ii^ — > e + ui + iO"*" , and then taking the 
limit (J — > 0, we obtain 



1 /■°° 

- / Auidu2A^{ui,U2) 

^ J —00 



X 



1 



+ 



1 



+ 



e + a;2 - iO+ € + 0Ji + iO+ 
da;ida;2 ^+(^1,(^2) 

1 1 



a;i — 6<;2 (e + W2 — iO^)(e + '^i + iO+) 



+ 



e — UJ2— iO+ e — 0J-1 + iO+ / a;i — u;2 



-/3a;i _ p-/3w2 



+ 



(e-u;2-i0+)(e-cc;i+i0+) 



(A-6) 



Note that A±{u2,uji) = A'^{L0i,iJ2) due to the time-reversal symmetry, and thus $1^1 (e,e) is real. 
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Fig. 3. Three analytic regions of the current vertex Aij;^^^^ [z, + w). 
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Fig. 5. Schematic picture of the model. 
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8/t 8/t 

Fig. 9. Transmission probability for N = 2. The temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 ■ ■ •, and 
(d) 0.5. The onsite interaction U/{2nt) is taken to be ( — ) 0.0, (-o-) 0.5, and (-•-) 1.0. The horizontal axis e is 
measured in units of t. Here T/t — 0.75. 
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£/t 8/t 

Fig. 10. Transmission probability for A'^ = 3. The temperature T /t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 ■ ■ ■, and 
(d) 0.5. The onsite interaction U/{2nt) is taken to be ( — ) 0.0, (-o-) 0.5, and (-•-) 1.0. The horizontal axis e is 
measured in units of t. Here F /t — 0.75. 
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£/t S/t 

Fig. 11. Transmission probability for = 4. The temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 ■ ■ •, and 
(d) 0.5. The onsite interaction U/{2nt) is taken to be ( — ) 0.0, (-o-) 0.5, and (-•-) 1.0. The horizontal axis e is 
measured in units of t. Here T/t — 0.75. 
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S/A E/A 

Fig. 12. The bubble (o) and vertex (•) contributions to the transmission probability (solid line) for TV = 1. The 
temperature T/A is taken to be (a) 0.0, (b) 0.1, (c) 0.333 • ■ •, and (d) 1.0. Here [//(ttA) = 2.0. 
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8/t S/t 

Fig. 13. The bubble (o) and vertex (•) contributions to the transmission probability (solid line) for TV = 2. The 
temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 • • •, and (d) 0.5. Here U/{2%t) = 1.0, and T/t = 0.75. 



24 




Fig. 14. The bubble (o) and vertex (•) contributions to the transmission probability (solid line) for TV = 3. The 
temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 • • •, and (d) 0.5. Here f//(27ri) = 1.0, and V/t = 0.75. 
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Fig. 15. The bubble (o) and vertex (•) contributions to the transmission probability (solid line) for TV = 4. The 
temperature T/t is taken to be (a) 0.0, (b) 0.2, (c) 0.333 • • •, and (d) 0.5. Here U/{2-Kt) = 1.0, and T/t = 0.75. 
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Fig. 16. Conductance vs temperature: U/{2nt) is taken to be (a) 0.5 and (b) 1.0 for N = 2 (□), A'^ = 3 (•), and 
N — 4 (O). The noninteracting results are also plotted for N = 2 (dotted line), = 3 (solid line), and = 4 
(dashed line). Here T/t = 0.75. 
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